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The effects of dust charge variations on low-frequency wave modulations in an electronegative dusty plasma 
are investigated. The dynamics of the modulated wave is governed by a nonlinear Schrodinger equation with a 
dissipative term. The dissipation arises due to the nonsteady (nonadiabatic) dust charge variations. Theoretical 
and numerical investigations predict the formation of dissipative bright (envelope) and dark solitons. The 
nonsteady charge-variation-induced dissipation reduces the modulational instability growth rate and introduces 
a characteristic time scale to observe bright solitons. Results are discussed in the context of electronegative dusty 
plasma experiments. 
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I. INTRODUCTION 

Electronegative dusty plasmas are very much common in 
modern technology (e.g., semiconductor material processing, 
development of nanomaterials, and for the efficient deposition 
of quality thin films) [1-6] and also are applicable in 
astrophysics (e.g., in the Earth's lower ionospheric plasmas 
and cometary plasmas) [7-9]. The dust grains produced by the 
chemical reaction in gas- or plasma- surface interaction in such 
plasmas can significantly affect the local as well as the global 
discharge characteristics [4,10]. 

In most of the astrophysical and laboratory plasmas, the 
electron number density decreases because of its attachment 
on the dust grain's surface. The presence of negative ions 
significantly affects the magnitudes of different plasma char- 
acteristic parameters as well as the dust charging process by 
modifying the surface potential acquired by a dust grain [11]. 
Actually, the thermal velocity of negative ions is lower than 
that of electrons, and as a result the presence of negative ions 
reduces the dust surface potential in an electronegative dusty 
plasma [12-14]. This causes the reduction of the magnitude of 
the dust charge. These modifications affect the properties of the 
low-frequency natural linear modes in a dusty plasma [15,16]. 
The nonlinear localized wave structures (e.g., a double layer) 
are also observed [13] in an electronegative dusty plasma. 

However, in a real dusty plasma the charge (qj) on dust 
grains play as an extra dynamical variable which can vary 
continuously with time due to continuous impact by electrons 
and ions (depending on the surrounding plasma environment). 
The instantaneous dust charge variation is known as steady- 
state (adiabatic) charge variation [17,18]. In this case, the 
dust charge instantaneously reaches its equilibrium value at 
each space-time point determined by the local electrostatic 
potential. However, the noninstantaneous charge variation 
causes an anomalous dissipation in a dusty plasma, which is 
known as nonsteady (nonadiabatic) charge variation [22] . This 
nonsteady charge variation causes a collisionless, non-Landau 
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damping of linear wave modes [15,19,20]. In the nonlinear 
regime, this anomalous dissipation leads to the formation of a 
shock wave in a dusty plasma [21-24]. In the long- wavelength 
limit, the shock is described by the Korteweg-de Vries Burgers 
equation [22,23]. The presence of negative ions modifies the 
nonlinear localized structures [25] and also reduces the shock 
strength by reducing the dust surface potential [26] . 

Moreover, in the short- wavelength limit, the low-frequency 
mode, e.g., a dust acoustic wave (DAW) [27], is modulated 
due to a small modulational perturbation in a dusty plasma 
[28-31]. A slow parallel modulation of a finite amplitude 
monochromatic plane wave can grow and in some limit leads 
to the formation of a bright (envelope) soliton. The physics 
behind this modulational instability is the interplay between 
nonlinearity and dispersion or diffraction. The nonlinearity 
mainly originates from the ponderomotive force. Some locally 
low density regions (cavity) are formed in a uniform dusty 
plasma due to this ponderomotive force. As a result, the 
DAWs are trapped within these cavities. This trapping of the 
plasma waves further enhances the ponderomotive force, and 
finally a growing envelope (bright) soliton is developed from 
a monochromatic plane wave in the unstable region. In the 
stable region, DAWs produce dark solitons or envelope holes. 

In this paper, we will investigate the DAW modulations 
in an electronegative dusty plasma in the presence of dust 
charge variations. In particular, we are interested to see how 
the negative ions and nonsteady dust charge variations affect 
the DAW modulations at short wavelengths. We consider the 
dust fluid equations to describe the model. We also adopt 
the well-known reductive perturbation technique (RPT) to 
study the DAW modulations. The dynamics of the modulated 
DAW is seen to be governed by a dissipative (here damped) 
nonlinear Schrodinger equation (NLSE). The damping arises 
due to the nonsteady dust charge variations. Both analytical 
and numerical solutions reveal that this damping causes the 
soliton (bright and dark) amplitude to decay with time. The 
presence of negative ions reduces (increases) the magnitude of 
the soliton amplitude (width). 

The paper is organized as follows. In Sec. II, we present 
the complete set of fluid equations to describe the model. The 
damped NLSE using RPT is derived in Sec. III. The effects 
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of dust charge variations on wave modulation are discussed in 
Sec. IV. The soliton solutions and the effects of dust charge 
variations on these solutions are discussed in Sec. V. The 
results are briefly discussed in Sec. VI. Section VII contains 
the concluding part. 

II. PHYSICAL ASSUMPTIONS AND BASIC EQUATIONS 

We consider an unbounded, unmagnetized, and collision- 
less electronegative dusty plasma. The plasma constituents 
are the electrons, singly charged ions (positive and nega- 
tive), and negatively charged dust grains of uniform size. 
The plasma is assumed to be in the undisturbed uniform 
state at — oo, where electrostatic potential = 0, electron 
number density n e = n e o, positive (negative) ion density n± = 
n±o, dust density rid = ndo, and dust charge q d = — Zd^, so 
that the quasineutrality condition n e o + Zd^do + ^-o = is 
satisfied. 

The charge on the dust grain varies continuously with 
space (x) and time (t). The temperature of the dust grain 
is very low compared to that of the electrons (T e ) and ions 
(7±, positive and negative). So, the dust grains are effectively 
cold with respect to the electrons and ions. The dust grains 
are moving with fluid velocity U. It is convenient to express 
all the variables in the nondimensional form before going 
to the details of the basic formalism of the problem. For 
this purpose, we introduce the following normalizations: 
O = ec/)/T e , N = rid/rido, N e = n e /n e o, N± = n±/n±o, x = 
x/k D , t = a) pd t, U = U/c da , and q d = q d /Zde = -1 + Qd , 
Qd is the fluctuating dust charge. Here, cD pd (=y/z 2 d e 2 n d o/€ m d ) 
is the dust plasma frequency, A (= ^/€oT e /yn+oe 2 ) is the 
plasma Debye length, c da (= y/zdT e /a d md) is me dust acous- 
tic speed, a d = z d n d o/yn + o, 8+ = n e0 /n +0 , 8- = n- /n +0 , 
o± = T±/T e , A = (1 — <5_ — y = 8 + + l/or+ + 
and z = Zd^ 2 /47t€ordT e is a dimensionless dusty plasma 
parameter (the ratio of the electrostatic energy of a dust grain 
of radius to the electron thermal energy). Hereafter, we will 
be using these new variables and removing all the bars for 
simplicity of notations. 

We are interested to study the low-frequency DAW mod- 
ulations, and therefore we neglect the inertia of the electrons 
and ions (positive and negative) compared to the dust grains. 
On this slow time scale, the electrons and both the ions are 
in local thermodynamic equilibrium, so that their densities are 
modeled by the following Boltzmann distributions: 

n e = n e oe®, n + = n + oe~® /a+ , ri- = ri-oe® /a ~ . (1) 

The normalized fluid equations for the dynamics of a DAW 
in a one-dimensional (ID) planar geometry in electronegative 
dusty plasma are 



dN d 

— + — (NU) = 0, 
at ax 

dU dU (Q d -l)d<& 

— + U — = -— — , 

dt dx ad dx 



(2) 
(3) 



and 



Y = + ^ ~ ~ AN(Qd ~ 1} ' W 



The other dynamical variable, the dust charge, is determined 
by the following normalized dust charging equation: 



oQd | v dQ d 
dt dx 



(l+OT + +/X 2 i )L 

— A+e^ +zQd) — A_e 



_ zQ d 
z + Cf + 



e a + 



(5) 



The charging frequency v ch is given by [26] : 



r d ^+(z + flr+)(l+or + + At 2 ) 
v ch = -^=T7 n > ( 6 ) 



2^ V t - 



ZCf+fr 



'ch 



where cD p+ (=^/n+oe 2 /€om+) is the positive ion plasma fre- 
quency, and V t+ (= ^T + /m + ) is the positive ion thermal 
velocity. The other coefficients A+, A_, fa, and /x^ are given 
in Appendix A [See Eqs. (A1)-(A3)]. 

III. NONLINEAR EVOLUTION EQUATIONS 

In order to study the propagation characteristics of a 
nonlinear DAW in the short-wavelength limit, we introduce 
the following slow space and time scales [32]: 



§ = €(x — v g t), 



X = € t, 



(7) 



where v g is the group velocity of the DAW and € is a small 
parameter that characterizes the strength of nonlinearity. This 
transformation separates the system into a slowly varying part 
associated with the amplitude of the wave and a rapidly varying 
part, which is the phase of the wave. The operators 3/3 1 and 
3/3x then take the forms 



and 



d 
dt 



3 

dx 
3 

> Yt 
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to account for the slow variations of DAW amplitude. 

The dependent variables are expanded in powers of € in the 
following way: 



"TV" 




"l" 


U 







Qd 







_ o _ 




_()_ 



ql n \i;,T) 



exp[i(kx — cot)l]. 



The reality condition for all variables is S 



in) 



c(») _ n ( n ) 



Lin) 



(8) 

(S\ n) f, where 



In modern technological experiments (e.g., Ar + — F ) 
and also in electronegative laboratory dusty plasmas (SiH^ — 
SiHp [2,3,10,15], the ratio co pd /v ch - O(10~ 2 - 1(T 3 ), 
which is small. These values of co pd /v C h permit us to consider 
the scaling 



0(€ 2 ), 



(9) 
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for the consistent perturbation. Using the stretching Eqs. (7) 
and (8) and scaling Eq. (9), we expand Eqs. (2)-(5) to obtain 
the nth-order reduced equations (B 1)-(B4) (details are given in 
Appendix B). For the first order (n = 1 and / = 1), we obtain 
the following relations: 



(1) 
con\ 


= ku[ l \ 


(10) 


(i) 
a d oou\ 


= -k^\ 


(11) 




= A(,<» -»?>), 


(12) 



and 



(13) 



qr = -Mi 

These four relations (10)— (13) self-consistently determine the 
dispersion relation, 

Ak 2 

(14) 



a d [y(l+k 2 ) + A& h ]' 

of the DAW in an electronegative dusty plasma. In the second 
order (n — 2 with / = 1), we actually get the corrections of 
the first-order quantities, which are written in the following 
matrix form: 



2iky iyco 



A A 
r k 2 k 



T 



(2) 



(15) 



where the superscript T represents the transpose of a matrix. 
The second order also provides the following compatibility 
condition: 



doj 



CO 

~k 



y + A# 



*ch 



y(l+k 2 ) + Ap ch ]' 



(16) 



In the above relations we use the conditions S$p =0 and = 
for |/| > 1. To the next order (n = 2, 1 = 2), we obtain 
the following second-harmonic modes (generated due to the 
self-interactions among modes): 



[0f,4 2) , 



■ (2 \q^] T = [A, B, C, 



(AhA + ft)] r 0{ 1)2 . 



(17) 



The expressions for A, ^2, B, and C are given in Appendix A 
[see Eqs. (A4)-(A7)]. The zeroth-harmonic modes that appear 
due to self-interaction among the modulated carrier wave 
are determined from the / = component of the third-order 
(n = 3) relations and are found to be given by 

' 1 ( 2k3 . 



2 1 

ot d oo 



(-/3 ch D + 2p 2 ) 



(18) 



The expressions of D and E are given in Appendix A [see 
Eqs. (A8) and (A9)]. Finally, we substitute the above derived 
relations into the third-order (n = 3) relations and obtain the 
following nonlinear Schrodinger equation (NLSE) for 0^ 
[= ir] : 



. d\j/ 



d 2 f 



(19) 



with a dissipative term T C h. In this NLSE, the group dispersion 
coefficient P [for the expression, see Eq. (All)] is related to 
the curvature of the dispersion relation co(k) [Eq. (14)], which 
is always negative for all wave numbers k. The coefficient of 
the nonlinear term Q [for the expression, see Eq. (A 12)] is 
related to the nonlinear frequency shift. The dissipative term 



r c h = 



2 12 



2c^[y(l+£ 2 ) + A& h ] 2 



^ ) (20) 



is related to the dust charge variations. This expression shows 
that the dissipative term T C h disappears for &> P d/v C h ~ 0, as 
it should because this implies that the left-hand side of the 
charging equation (5) equals zero, i.e., I e + /+ + /_ = 0. In 
this case the dust charge Q d is in the local equilibrium state, 
providing an approximate description of the charge state of a 
dust grain, and produces no dissipative effect (adiabatic dust 
charge variations). Also note that r ch = for A(=l — 8- — 
8+ = z d n d o/n+o) = 0; this is expected as A = implies the 
absence of charged dust grains and consequently no charge 
variations. Hence, the dissipative term T C h arises only due to 
the nonsteady (nonadiabatic) dust charge variations. 

IV. EFFECTS OF DUST CHARGE VARIATIONS 
ON MODULATIONAL INSTABILITY 

At very low frequencies, the DAW is nearly at zero 
frequency and Eq. (19) can be analyzed for its stability. In the 
above damped NLSE [Eq. (19)], the term —Q\if\ 2 plays the 
role of a potential energy. The local maxima of this potential 
energy acts as an effective potential well in the plasma. The 
physics of this phenomena can be explained in the following 
manner. The high-frequency DAW reflects from the regions 
of high density, and the amplitude-dependent ponderomotive 
force forms a cavity (low-density regions) in the plasma. The 
DAW becomes trapped within these density-depleted regions 
(cavity) made by the wave itself. If a wave of amplitude ifr is 
trapped in a well, its energy | \J/ \ 2 will concentrate at the bottom 
of the well, but the energy concentration makes the well deeper 
by making this energy even larger. Thus the formation of the 
density depletion is an unstable physical process that occurs 
due to the energy localization in the medium. This process 
leads to an instability known as modulational instability (or 
caviton instability in the case of one dimension) [33,34]. Here, 
we briefly discuss this instability for DAW in electronegative 
dusty plasma and the effects of charge variations on this 
instability. 

Let us assume that, in the presence of nonadiabatic dust- 
charge-variation-induced dissipation, the amplitude evolution 
equation (19) possesses the following plane- wave solution: 



f = VoMexp 



8(x')dx' 



(21) 



where ^o(^) and 8(r) are the amplitude of the pump carrier 
wave and the nonlinear frequency shift in the presence of 
dissipation. Substitution of this plane- wave solution in Eq. (19) 
yields the following two equations: 



dr 



+ r ch ^o = f Q {x) = ^00 exp(-r ch r), (22) 
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and 



8(r) = -Q\Mr)\ 2 



-£IVoo| 2 exp(-2r ch r), (23) 



where ^oo is a real constant. Also note that VoW ^ as 
r — >► oo, which implies that x/tq is stable. For stability analysis, 
we consider the perturbation about this stable solution in the 
following standard procedure: 



f = [Vo(t) + f(^r)] exp -i / 8{T')dx' 



f 

Jo 



(24) 



where \jr{\ \jr V^o) is the perturbed amplitude of the modu- 
lated wave. 

After the substitution of Eq. (24) in Eq. (19), one can show 
that the perturbation satisfies the following linearized equation: 

+ P S + 0l^o| 2 OA + #*) + = 0, (25) 

where tfr* is complex conjugate of yjr. Let us put \[r = fy R + 
/^/, and then the above equation satisfies the following two 
coupled equations: 



and 



3r 



a§ 2 



(26) 



(27) 



Finally, the space-time dependence of the perturbation of 
the form \j/ ~ exp(/#) [where ft = k% — J* a)(z)df is the 
modulated phase with k(<^ k) and <£>(^C &>) being the mod- 
ulation frequency and wave number, respectively) yields the 
dispersion relation: 

(& + zT ch ) 2 = P 2 k 4 - 2PQ\xlf \ 2 P. 

This dispersion relation shows that the nonsteady (nonadia- 
batic) dust-charge- variation-induced dissipation provides the 
usual damping of the wave. The system is stable for P Q < 0. 
However, there is a possibility of instability if P Q > 0: both P 
and Q are of the same sign, here P < for all k [Eq. (All)], 
so to make PQ > we must have Q < 0. In this case of 
nonadiabatic dust charge variations, the instability could occur 
if 

k 2 < k 2 cr = (^pj | Vol 2 = (jP) |^oo| 2 exp(-2r ch r), 
provided with these values of k the inequality 



must hold. This determines the maximum time, 

r m;iX = nz In i 



1 . f2PQk 2 \xlfr a 



2r, 



ch 



^00 1 

r c 2 h + p^4 



(28) 



to observe instability. Thus the instability growth will cease 
for r ^ r max . For PQ > 0, by setting 6b = zT, we obtain the 
dispersion relation, 

(T + T ch f = 2PQ\^\ 2 k 2 -P 2 k\ 



for the instability growth rate. The maximum growth rate 
is obtained by taking the derivative of the above dispersion 
relation with respect to P, and setting this to zero. This yields 



f ) m 2 



Q 



| ^oo I exp(-2r ch r), 



and, with these values of £^ ax , me maximum growth rate 
becomes 

r max (nonadiabatic) = V2| g| |^o| 2 - r ch 

= V2|G||Voo| 2 exp(-2r ch r) - r ch . 

(29) 

The configuration is unstable when a/2 I Q 1 1 Vo 1 2 > T c h. 

In the case of adiabatic charge variations, &> P d/v C h = 0, i.e., 
T ch = 0, and there is no dissipation in the system. In this case, 
one can easily find the dispersion relation, 

cb 2 = P 2 k 4 -2PQ\xlf 00 \ 2 k 2 , 

of the modulated wave. The instability criteria for PQ > 
of the modulated wave then becomes k 2 < k 2 r = (2Q/P) 
| Vool 2 - In this case, 



Fmax (adiabatic) = lim r max (nonadiabatic) 

r ch ^o 

= V2\Q\\xlf 00 \ 2 



(30) 



is the maximum growth rate. Comparing Eq. (30) with that of 
the nonadiabatic case [Eq. (29)], we obtain 

I^max (adiabatic) > r max (nonadiabatic), 

as T ch > 0. 

All the above investigations reveal that the nonadiabatic 
dust-charge- variation-induced dissipation reduces the insta- 
bility growth and therefore introduces a characteristic time 
scale to make the instability occur. Thus the nonlinear DAW 
becomes modulationally unstable within the time interval 
[0, r max ) for P Q > 0. The nonlinear saturation of such DAWs 
can give rise to bright (envelope) soliton structures. Therefore, 
the value of r max provides an estimation of a characteristic time 
to observe a bright soliton in a dusty plasma with nonadiabatic 
dust charge variations. However, in the adiabatic limit (r c h —> 
0) r max —> oo. Thus, in the case of steady-state (adiabatic) 
charge variations, there does not exist any characteristic time 
scale to observe the bright soliton. In the stable region P Q < 0, 
it is well known that a finite amplitude DAW can not form a 
envelope soliton [35]; instead, it forms a dark soliton structure. 

We now investigate the group dispersion coefficient P and 
nonlinear coefficient Q and also the instability conditions 
numerically for wave number k. The computations are car- 
ried out for the representative parameters of the Ar + — F~ 
electronegative dusty plasmas, used for ultrafine and highly 
selective etching of polysilicons [2,10,15]. The parameters 
used in the computations are summarized as follows: n+o ~ 
4 x 10 17 m- 3 , n d0 ~ 4 x 10 13 m- 3 , T e = 2eV, T_ = O.leV, 
T + = 0.2eV, m_ - 3.2 x 10~ 26 kg, m+ - 6.69 x 10~ 26 kg, 
and the dust mass density ~ 1.5 x 10 3 kgm -3 . The dusty 
plasma parameter z is computed for two different dust 
grain radii and four cases of negative ion content <$_(= 
n-o/n+o)- The cases are [15] (a) rj = 0.1 /xm, 5_ =0.1; 
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FIG. 1. (Color online) Variations of Q with wave number k for 
different plasma parameters. The plasma parameters for different 
figures are (a) r d = 0.1 fim, 8- = 0.1, z = 3.3, and 8 + = 0.854; (b) 
r d = 0. 1 fi m, 8_ = 0.4, z = 2.98, and 8+ = 0.56; (c) r d = 1 /i m, 
8_ = 0.1, z = 2.89, and 8+ = 0.498; and (d) r d = 1 /i m, <3_ = 0.4, 
z = 2.42, and 8+ = 0.337. 

(b) r d = 0.1 /xm, 5_ = 0.4; (c) r d = 1 /xm, 8 = 0.1; and (d) 
r d = 1 /im, 5_ = 0.4, respectively. 

Different curves (a), (b), (c), and (d) in Fig. 1 show that the 
magnitude of the coefficient of nonlinearity Q increases with 
the increase of negative ion content. However, a comparative 
study between the different curves shows that Q decreases for 
larger dust grains. The magnitude of the coefficient of group 
dispersion P increases with the increase of the proportion 
of the negative ions for larger dust grains. But comparative 
study shows that the difference is insignificant for smaller dust 
grains. These are shown graphically in Fig. 2. 

The nonlinear DAW becomes modulationally unstable only 
if Q < (as P is always negative). Figure 1 reveals that 
Q < for the wave number k < 0.5 and also for k > 1. 
Thus DAW becomes modulationally unstable in the very- 
small-wavelength perturbation and also in the sufficiently- 
long- wavelength perturbation. The nonsteady dust-charge- 
variation-induced dissipation T C h decreases with the increase 
of the negative ion concentration, as shown in Fig. 3. This 
figure also shows that the dissipation is very low for larger 
dust grains. 

V. SOLITON SOLUTIONS 

In the above damped NLSE [Eq. (19)], the group dispersion 
coefficient P [Eq. (All)] is always negative, but the nonlinear 
coefficient Q [Eq. (A 12)] can be both positive and negative 
depending on the wave number k. Thus, to find the soliton 
solutions of the damped NLSE [Eq. (19)], we recast the 
equation in the following normal form: 



l ~di 



l aV 

2~W 



(31) 



where f = | Q|r,f = § VI Q\JW\ and #) = r^/|fi|. 
This equation exhibits a bright or envelope soliton for k = — 1 , 
whereas it exhibits a dark soliton for k = +1. Equation (31) is 
solved numerically for k = ±1, and the solutions are shown 
graphically in Figs. 4 and 5. However, in the proceeding 
subsections, we have solved this equation analytically. 
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FIG. 2. (Color online) Variations of group dispersion coefficient 
P with wave number k for different plasma parameters. The plasma 
parameters for the different figures are the same as in Fig. 1 . 

A. Bright soliton: adiabatic charge variation 

In the case of adiabatic dust charge variations, co p d/ v ch ~ 
0, and then, in Eq. (31), we can put e(^) = 0. Thus, in the 
adiabatic limit, we have the usual NLSE with k = — 1 . Let us 
assume a solution of the form \//(^,f) = p(f,f ) exp[/<p(f,f )] 
and then solve the ordinary differential equations for cp and 
p subject to the boundary condition p^0as§^±oo. We 
finally obtain the following single envelope (bright) soliton, 
which in terms of actual parameters reads as 



^r(£,r) = a cosh 



x exp 



Vi\PQ\xr) 




(a z -X l )r 



(32) 



where a and / are two soliton parameters in which a is the 
amplitude of the soliton. This shows that the disturbances 
resemble the soliton shape with an exponential factor making 
oscillation between a maxima and a minima. The resultant 
structure is the envelope excitation of DAW. The graphical 
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FIG. 3. (Color online) Variations of nonsteady dust-charge- 
variation-induced dissipation r ch with wave number k for different 
plasma parameters. The plasma parameters for the different figures 
are the same as in Fig. 1. 
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FIG. 4. (Color online) Bright soliton in the presence of dissipation 
due to nonsteady dust charge variation. Numerical solution of Eq. (31) 
for ic = -1, where f = \ Q\x andf = \~J\Q\p\F\. 



representation of this single bright soliton is shown in Fig. 6 
for different plasma parameters. 



B. Bright soliton: nonadiabatic charge variation 

In the case of nonadiabatic dust charge variations, 
&>pd/^ch 7^ 0, and if this dissipation is weak [here it is 
indeed weak, as &> P d/v C h ~ 0(e 2 )], we can solve the above 
damped NLSE [Eq. (19)] perturbatively by taking s{\jf) as 
a small perturbation. Therefore, to solve it, let us take the 
general solution of the perturbed soliton in the following form 
[36,37]: 



^(§,r) = a(r)cosh + b(r)]} exp[/§x( r ) 



icr(r)], 
(33) 



where a, b, cr, and / are the soliton parameters. Finally, 
applying the conservation laws for the NLSE (conserved 
integral relations [36,37]), we obtain the following bright 
soliton (envelope soliton) in the presence of nonsteady dust 
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FIG. 5. (Color online) Dark soliton in the presence of dissipation 
due to nonsteady dust charge variation. Numerical solution of Eq. (31) 
for a: = +1. The independent variables are the same as in Fig. 4. 
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FIG. 6. (Color online) Bright (envelope) soliton. Numerical 
solution of Eq. (19) with r ch = for PQ > with k = 0.2. The 
plasma parameters for the different figures are the same as in Fig. 1. 

charge variations, which in terms of actual variable reads as 
^r(£,r) = «oexp(— 2r c hT)cosh _1 



a exp(-2r ch r) / +^2\P Q\ Xo t) 



x exp 




(34) 



where ao and xo are the initial values of a and x » respectively. 
In the limit T C h — ► 0, we recover the previous result [Eq. (32)] 
(with ao = a and xo = X)- It is clear that as time elapses 
the amplitude of the bright (envelope) soliton decreases 
exponentially with decay rate ~ 2r ch . The numerical solution 
in Fig. 4 also shows a similar nature. Thus the nonadiabatic 
dust charge variation has a damping effect on the bright soliton 
structure in a dusty plasma. 

C. Dark soliton: adiabatic charge variation 

As before, in this case we have the usual NLSE with 
k = + 1 . To obtain the dark soliton solution, proceeding as be- 
fore, subject to the boundary conditions p — >► uq (a constant), 
we obtain the following single dark soliton [36,37,39]: 



W,t) = u [rj tanh(0 + iy/ 1 - rj 2 ] exp (i | Q\u 2 t) , (35) 



where f = w ^(v1 Q\/2\P\$ + u z \ Qy 1 - r] 2 ). The param- 
eters uo and r] represent the amplitude and depth of the dip of 
the dark soliton. The single dark soliton is shown graphically in 
Fig. 7 for different plasma parameters. Dark solitons are called 
gray solitons when \ rj\ < 1 . On the other hand, the | rj \ = 1 case 
corresponds to the black solitons or stationary shock solution 
given by 

\K$ 9 r) = uq tanh[ Mo V I Q\/2\P\^]Gxp(i\Q\ulr). (36) 

In this case an input pulse exhibits an intensity "hole" at the 
center, and thus, sometimes, it is called an envelope hole. 
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FIG. 7. (Color online) Dark soliton. Numerical solution of 
Eq. (19) with r ch = for PQ < with k = 0.5. The plasma 
parameters for the different figures are the same as in Fig. 1. 

D. Dark soliton: nonadiabatic charge variation 

In the presence of weak dissipation s{yjr), we assume 
\j/(^,z) = ^[wo(t),77(t);£,t] using the above solution (36) 
and then applying the conservation laws for the NLSE 
(conserved integral relations [38,39]), we obtain the following 
dark soliton: 



1r($,r) = w (r)Mr)tanh[f ( r )] + iy/l - t](t) 2 
x exp(/|g|wo(r) 2 r), 



(37) 



in the the presence of nonsteady dust charge variations. 
In the above f(r) = u (t)t](t)(^\Q\/2\P^ + u (r)r\Q\ 



r](r) 2 ) and 
u (r) = w(0)exp(-r ch r), 



^(r) = ^(0)exp(-r ch r), (38) 



where u(0) and r](0) are the initial values of uq and respec- 
tively. The dark soliton amplitude also decays exponentially 
with time with decay rate ~ T C h due to the nonsteady dust 
charge variations. The numerical solution in Fig. 5 also shows 
a similar nature. In this case the decay rate is half of that of 
the bright soliton. In the limit T C h — > 0, we recover the usual 
single dark soliton [Eq. (35)]. 

VI. DISCUSSION 

In the preceding sections, we see that the dynamics of 
the modulated DAW in charge varying electronegative dusty 
plasma is governed by a damped NLSE. The damping term 
appears due to the nonsteady (nonadiabatic) dust charge 
variations. To study the effects of nonsteady dust charge 
variations on DAW modulations, we have assumed that 
&>pd/^ch ~ 0(€ 2 ). This assumption is valid if the dust oscil- 
lation period is comparable with the dust charging frequency, 
i.e., the conditions &> P d/v C h <3C 1, ^ are met. The specified 
Ar + — F~ laboratory plasma parameters yield co V d/ v ch ~ 
1.4 x 10 -2 for smaller dust grains (rj =0.1 fim). Thus, the 
condition for nonadiabaticity is indeed satisfied. However, for 
larger (r^ > 1 /xm) dust grains the ratio becomes <4 x 10 -4 
and the conditions could be violated. In this case, one can 
consider the adiabatic dust charge variations. 

The present investigation shows that the DAWs are trapped 
in a hole created by the wave itself in electronegative dusty 



plasma and they become modulationally unstable in the 
region PQ > 0. In this unstable region DAWs admit a bright 
(envelope) soliton, as shown in Fig. 6. In the stable region 
P Q < 0, DAWs admit a dark soliton, as shown in Fig. 7. The 
presence of negative ions decreases the magnitude of the bright 
soliton amplitude (depth of the dip of the dark soliton), and, 
as a consequence, the width of both solitons increases. Also, 
with the increase of the dust grain radius, the magnitude of the 
soliton amplitude (depth of the dip) increases, whereas width 
increases [see the comparative study between the solitons (a) 
and (b) and the solitons (c) and (d) in Figs. 6 and 7]. 

In the case of steady-state (adiabatic) dust charge variations, 
the dust charge instantaneously reaches its equilibrium value 
and thus does not produce any dissipative effect. However, 
in the case of nonsteady dust charge variations the growth 
rate of modulational instability is reduced. Actually, due to the 
delayed charging (nonadiabatic dust charge variations), i.e., the 
nonzero value of &> P d/v C h, there is an extra contribution from 
the dust charge Qd through the charging equation (5). This 
causes a decay of the dust charge magnitude and, consequently, 
plays the role of the dissipative mechanism on the dust fluid, 
causing dust motion by decreasing the driving force (oc Qd)- 
This dissipation produces the damping effects on collective 
phenomena in a dusty plasma. The damping effect opposes the 
energy localization in the medium (that causes the instability), 
and as a result the instability growth rate is reduced. 

The nonsteady dust charge variations also introduce a 
characteristic time r max to observe modulational instability, 
i.e., to observe envelope solitons, which cannot be found in 
dusty plasma with a constant dust charge [29], and also for 
steady- state charge variations. This is a unique feature of a 
truly dusty plasma. In the presence of dissipation, the analytical 
solutions show that the nonlinear DAW can propagate in the 
form of bright and dark solitons with decaying amplitude (and 
amplifying width). The decay rate of the bright soliton is twice 
the decay rate of the dark soliton. The numerical solutions also 
reveal the same nature as that shown in Figs. 4 and 5 and thus 
confirm the analytical results. 

It should be noted that the physical parameters for the 
adiabatic assumption are different from those of the nona- 
diabatic assumption. As a result, both bright (envelope) and 
dark solitons for the adiabatic case are quantitatively different 
from those of the nonadiabatic case. 

VII. CONCLUSION 

On the basis of dust fluid equations, we have analyzed the 
effects of dust charge variations and negative ions on DAW 
modulations in a dusty plasma. Our investigation shows that 
there is a possibility of the trapping of DAWs in a hole created 
by the wave itself in the medium. The parameter region, where 
DAWs become modulationally unstable, is specified. In the 
presence of nonsteady (nonadiabatic) dust charge variations, 
the dynamics of the modulated wave is governed by a damped 
NLSE. The modulated wave exhibits both dissipative (here 
damped) bright (envelope) and dark solitons. Negative ion 
content reduces (increases) the soliton amplitude (width). It 
should be noted that the balance between dispersion and 
nonlinearity gives rise to a similar type of bright (dark) 
solitons described by the NLSE in optical fiber communication 
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systems, provided that the pulse propagates in the anomalous 
group-velocity dispersion regime. The dissipative bright (dark) 
solitons are also observed in optical fiber, where the dissipation 
occurs due to the fiber losses in the system [36-39]. Here, the 
decay rate of the dark soliton is half of that of the bright soliton, 
which is similar to the observed dissipative optical soliton in 
the communication system [36-39]. 

Finally, we stress that the results of our investigation could 
be useful in understanding the features of modulated DAW 
packets in a dissipative electronegative dusty plasma. Recently, 
a dissipative dark soliton was observed in a strongly coupled 
complex plasma [40,41]. In that experiment, the dissipation 
(damping) occurs due to the neutral drag, whereas here the 
damping is due to nonadiabatic charge variations. In this 
connection, we must mention that we have not encountered 
any experimental observations of bright (envelope) and dark 
solitons in a dissipative electronegative dusty plasma. It would 
be very interesting to look at this situation in a laboratory. We 
hope that, in the future, such a type of dissipative soliton could 
be observed in electronegative dusty plasma experiments. 
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APPENDIX B: OUTLINES OF CALCULATIONS 



Substitution of Eqs. (7)-(9) in Eqs. (2)-(5) yields the following nth-order reduced set of equations: 
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